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Abstract 

Years ago S. Weinberg suggested the " Quasi-Particle" method (Q-P) for iteratively solving an 
integral equation, based on an expansion in terms of sturmian functions that are eigenfunctions of 
the integral kernel. An improvement of this method is presented that does not require knowledge of 
such sturmian functions, but uses simpler auxiliary sturmian functions instead. This improved Q-P 
method solves the integral equation iterated to second order so as to accelerate the convergence of 
the iterations. Numerical examples are given for the solution of the Lippmann-Schwinger integral 
equation for the scattering of a particle from a potential with a repulsive core. An accuracy of 
1 : 10^ is achieved after 14 iterations, and 1 : 10^'' after 20 iterations. The calculations are carried 
out in configuration space for positive energies with an accuracy of 1 : 10^^ by using a spectral 
expansion method in terms of Chebyshev polynomials. The method can be extended to general 
integral kernels, and also to solving a Schrodinger equation with Coulomb or non-local potentials. 
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I. INTRODUCTION 

Negative energy sturmian functions (S-F), introduced in the 60'ties by M. Rotenberg 
found many useful applications, such as in the calculation of electron induced ionization 
collisions in the identification of resonances in nucleon- nucleus scattering {4], in the 

solution of a Schrodinger equation with non-local potentials js], for a separable representation 
of scattering t— matrices |6| and in the solution of three-body Faddeev equations [7]. For 
the applications involving long range Coulomb forces, the analytical expressions for the 
Coulomb Green's functions, initially developed by L. Hosteler and R. Pratt [s^ have found 
many applications [9]. sturmian functions form a complete, discrete set of eigensolutions 
of a Sturm-Liouville differential (or integral) equation, and form a complete set of basis 
functions. However the expansion of a wave function into a set of sturmian functions in 
ma„v cases does not „ we.. Q. and met..ods to improve t..e eonvergenee. sueh a. 
Fade approximations, have frequently been utilized [2[. The atom-atom potentials that 
occur in atomic physics calculation often have strong repulsive cores, in addition to a long 
range attractive part, and for such potentials it has been difficult to obtain a reliable set of 
sturmian functions. A method based on spectral expansions in terms of Chebyshev functions 
has been recently developed 12| that overcomes this difficulty, and hence the use of sturmian 
functions is again of interest. Nevertheless, for numerical calculations the expansions have 
to be truncated, and because of the slow convergence of sturmian expansions, an iterative 
method to correct for the truncation error becomes desirable. It is the purpose of this study 
to introduce such a method, based on the original Quasi-Particle method of Weinberg 13| . 
and, for the purpose of testing it, apply it to the scattering of a particle from a potential as 
described by the Lippmann-Schwinger (L-S) integral equation. Since, the method described 
here should be applicable to more general kernels of an integral equation, the present study 
is done in anticipation of solving the more complicated two-dimensional integral equations 
that occur in the solution of three-body equations in configuration space 10|], 1^. The 
expected accuracy is to be better than 6 to 7 significant figures, desirable for doing atomic 
physics calculations. By comparison, the solution of three-body equations for nuclear physics 
applications, done commonly in momentum space jl5|, achieve an accuracy not better than 
four significant figures [l6|. 

The present paper examines the convergence of several iteration methods for the solution 
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of the L-S equation describing the scattering of a particle with positive energy E by a 
local potential V in one dimension. The insights gained from this example is to guide one 
'or solving a more general integral equation. The method initially suggested by Weinberg 
I3I . also called the " Quasi-Particle" method Q-P, requires sturmian functions that are 
the eigenstates for the exact integral operator. The present method, denoted as method S2, 
differs from the Q-P method in that it uses auxiliary sturmian functions that are based on an 
auxiliary potential V that is not the same as the potential V. Since it may be computationally 
easier to obtain the auxiliary sturmians rather than the sturmians for the original integral 
operator method ^2 may prove advantageous for the case of very complicated integral kernels. 
For the L — S example, the auxiliary sturmians are used to approximate the scattering 
potential V{r) by a non-local separable representation Vi\i{r,r') in terms of auxiliary 
sturmian functions and subsequently correcting for the truncation error Vr = V — Vn by 
means of iterations. We find that the rate of convergence of the iterations does not depend 
significantly on the choice of V ^ provided that the range of V is sufficiently larger than the 
range of V . The sturmian functions in all these cases are calculated in configuration space 
for a positive energy so that the asymptotic form of the approximated wave function 
is the same as the exact wave function. When solving for a bound-state wave function, it 
is customary the represent the potential by negative energy sturmians that are real, but 
for scattering cases, positive energy sturmians should be preferable. An advantage of the 
sturmian expansion method over the Fourier-Grid method 17|] for positive energies is that 
the sturmian method emphasizes only the spatial region where the potential is non-negligible, 
the asymptotic part of the wave function being already incorporated into the sturmian basis, 
while in a Fourier-Grid method, the asymptotic part has to be obtained explicitly. 

Since an accuracy of the scattering solution to 1 : 10^ is envisaged, the sturmian functions 
are also required to have the same, or better, accuracy. Our approach for solving the Sturm- 
Liouville eigenvalue equation 121] utilizes a "spectral" expansion method into Chebyshev 
polynomials [l^, 13], that permits one to prescribe a desired accuracy. With this method 
one can obtain sturmian eigenvalues and eigenfunctions with an accuracy of, for example, 
1 : 10^^, and because of its stability, one can incorporate the effect of long-range tails of 
potentials, as was demonstrated in the calculation of the bound state eigenvalue of a Helium- 
H... di.er H. The present ca.cu.at,o,. a.e performed in c„nfi,n.ation space, and can 
be generalized to include asymptotic Coulomb functions. Contrary to what is often done at 
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negative energies, our approach does not expand the Green's functions into a separable set 
of sturmians, but rather expands the potential V{r) into such a set. The present method to 



2l| , in that it does 



calculate sturmian functions is similar to a method introduced previously 
not use a square well potential sturmian basis set in terms of which the desired sturmians 
are expanded, and, being based on a spectral method, is considerably more precise. This 
additional feature permits a more accurate study of the iteration convergence properties. 

In section 11 we present the formalism that defines the sturmian functions; section 111 
describes the approximate Q-P methods Si and S2 that extend Weinberg's Q-P original 
iterative method so that the use of sturmian functions for the integral kernel are avoided, 
and auxiliary sturmians are used instead; in section IV numerical examples are presented for 
the case that the integral kernel is of the Lippmann-Schwinger type, and section V contains 
the summary and conclusions. Appendices A, B, C and D contain further details. 

II. NOTATION AND EQUATIONS 

The one- dimensional Schrodinger equation for the L = radial wave function ip{r) is 

V^p = V^|J (1) 

where 

is the differential operator, E is the energy assumed positive here, V is the scattering po- 
tential. These quantities are in units of inverse length squared, and were transformed from 
their energy units by multiplication with the well known factor 2fi/h?. The wave number k is 
related to the energy hj E = k"^. The corresponding Lippmann-Schwinger integral equation 
for ip is 

POO 

^{r) = F{r)+ goir,r')V{r')^{r')dr'. (3) 
Jo 

The undistorted Green's function is given by 



6^o(r,r') = -^F(r<)xif(r>) (4) 

where (r, r/) = (r<,r>) if r < r' and (r, r/) = (r>,r<) if r > r'. In the present zero angular 
momentum case 

F(r) = sin(/cr); H{r) = cos(/cr) + ism{kr). (5) 



According to Eqs. ([3]) and (jl]) the asymptotic form of ip is 

^(r ^ oo) = F{r) + S H{r) 



(6) 



S = -\ I F{r')V{r')i){r') dr' . (7) 



with 

Near the origin '?/'(r — 0) = 0, because both F as well as the integral term in Eq. ([3]) go to 
zero. 

The sturmian functions s = 1, 2, ... obey Eq. ([T]) with V replaced by A 

I? = AsV (8) 

where V{r) is the sturmian potential chosen conveniently, that need not be equal to the 
scattering potential V{r), and the A's are the eigenvalues to be determined. The sturmian 
functions $s obey the boundary conditions 

$,(r->0)=0; <i>s{r ^oo) = is H{r), (9) 

where the constant is determined by the normalization of the sturmian function. The 
functions $s obey the integral equation version of Eq. ([8]) 

POD 

Vs^sir)= / goir,r')Vir')^sir')dr', s = 1,2,3,.... (10) 
Jo 

with the Green's function given by Eq. (jl]) for positive energies. The rj^s are eigenvalues of 
the operator QqV, they form an infinite set with point of convergence at 0, and are related 
to the A's according to ris = l/Ag. With the choice (jlj) and of the Green's function the 
sturmians $s obey the boundary conditions ([9]), and 



/ F(r)V(r)<!>sir)dr (11) 

■ Vs Jo 



k i]s Jo 

The $'s are not square integrable for positive energies or orthogonal to each other. How- 
ever, they are orthogonal if one includes the potential as a weight function 

POO 

<^s\V^s'>= / ^s{r) V ^s'{r)dr = for 'r]s r]s/. (12) 
Jo 

In the bra-ket notation above, < $s| is not the complex conjugate of as is usually implied. 
The validity of (|T2|) can be shown by replacing $5/) by As'QqV^s') in the integral ([T2|) . and 



subsequently replacing {^sV Qo by r]s{^s- The result is < '^'s' > = Vs-^s' < ^s\V ^s' >, 
and if rjsAs' ^ 1, the above identity becomes absurd unless < ^s\V ^s' > =0. Because of 
the completeness of the sturmian functions, one has the identity 

which can also be written as 

oo _ 

S(r - = E ^(O'^s(r) ^^^-^^^ (14) 

depending on whether the delta function applies to the left or to the right of the integrand 
in an integral. 

One can understand intuitively the properties of the $'s as follows 211]. As the index s 
increases, the corresponding values of increase, and hence the potential A^V^ increases in 
magnitude. If V is real and attractive and the real part of A^ is positive, then the real part 
of AgV becomes more attractive, and the corresponding eigenfunction $5 becomes more 
oscillatory inside of the attractive region of the well. So, from one s to the subsequent 
one s + 1 the eigenfunction acquires one more node inside of the well. According to flux 
considerations the imaginary part of A^V^ has to be positive, i.e., the well has to be emissive 



2l| . Near the origin the flux is ~ since $ ~ however asymptotically the outgoing form 



of the wave function produces a positive outgoing flux. This outgoing flux is generated by 
an emissive imaginary potential, exactly the opposite of the case of an optical potential, 
that absorbs flux. These properties will be verified in the numerical section below. 



A. Numerical properties of sturmian functions. 

In the examples given below, the scattering potential V is of the Morse type with a 
repulsive core near the origin (Vp), and the sturmian potential V is either also of the Morse 
type with the repulsive core suppressed (V5), but of the same range as the scattering potential 
Vp, or it is a potential of the Woods-Saxon type (Vws) with a range larger than that of 
either Vp or Vs. The Morse potentials are formed by the combination of two exponentials 

V{r) = Vm exp(-r/a + /3/a) x [exp(-r/a + (3/a) - 2] , (15) 

with parameters given in Table |T] , and the Woods-Saxon potential is given by 
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TABLE I: Parameters for two Morse- type potentials P and S. 
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FIG. 1: Three sturmian potentials as a function of radial distance, given by Eqs. (I15p and ()16p . 



Vws = Vo/{l-exp[{r-R)/a]} 



(16) 



with Vq = —5 fm~'^, R = 15 fm and a = 0.5 /m.The resulting dependence of these 
potentials on r is illustrated in Fig. [H 

The spectrum of the eigenvalues for the two Morse-type potentials, with k = 0.5/m~^ 
is illustrated in Figs. [2] and El Since the potential Vs is entirely attractive, the real parts 
of As are positive, while the imaginary parts are negative, in accordance with the argument 
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FIG. 2: The spectrum of the sturmian eigen values A for the Morse Potential Vs defined in Table 
Uand illustrated in Fig. [TJ The wave number is A; = 0.5 fm"^. 
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FIG. 3: Same as Fig. [2] for the potential Vp defined in Table HI This potential has a repulsive 
core. 




r(fm) 



FIG. 4: Sturmian eigenfunctions for the potential Vs defined in Table H] for a wave number k = 
0.5 fm^^. They are normalized such that asymptotically they all approach the outgoing Hankel 
function H{r) = cos{kr) + i sin[kr). 



given above. These sturmian eigenfunctions are obtained by an adaptation of Hartree's 
iterative method designed to obtain energy eigenvalues to the Schrodinger equation, and 
described in Appendix C. A hst of these eigenvalues precise to 11 significant figures is given 
in the Appendix D. Some of the eigenfunctions are illustrated in Figs. Il]and|5l that show 
that at large distances where the potential becomes small, these functions become linearly 
dependent. The $s used in the present examples are normalized such that asymptotically 
they equal the function H{r), with unit coefficient. 

Because potential Vp has both a repulsive and an attractive part, the eigenvalues fall 
into two categories. In category I the eigenvalues A have a positive real part and a negative 
imaginary part, and the corresponding eigenfunctions are large mainly in the attractive 



8 



_1l . . . 1 

5 10 15 20 

r(fm) 



FIG. 5: Same as Fig [5] for the imaginary part of the sturmian functions. 
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FIG. 6: Real parts of sturmian functions <I>s for the potential Vp, for k = 0.5 fm . This potential, 
defined in Table HI has a repulsive core. 

regions of the potential well. Examples are given in Figs. [6] and [71 In category II the real 
parts of A are negative so as to turn the repulsive piece of the potential near the origin into 
an attractive well, and the formerly attractive valley into a repulsive barrier. Examples of 
the corresponding sturmian for indices s = 4, 7, 10, 15, and two of these functions are shown 
in Figs. M and |9l For < r < 4 fm the function $7 has a magnitude less than 10~^. 
The sturmian for s = 10 is similar to that for s = 4 in that it is also large near the origin 
(with an amplitude of ~ 10^) and has a node near r = 1. The functions for s = 4 and 10 
are "resonant" in the radial region r G [0,4], while the one for s = 7 is non-resonant. For a 
larger energy the effect of the barrier for the functions of class II decreases, and the absolute 
value of the eigenvalues rjs = l/Ag decreases for s < 10, as is illustrated if Fig. [TOl As will 
be shown in section III the functions for which rjs < 1 play an important role in the Q-P 
method for the iterative correction of the truncation errors. Further, in the expansion of a 




FIG. 7: Same as Fig. [6] for the imaginary parts of the sturmian functions 
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FIG. 8: Sturmian function in category II for s = 4, for potential Vp and k = 0.5. 



wave function in terms of sturmians, such as described in section IV-A below, a dominator 
(1 —77s) is likely to appear. For the values of s for which real(?7s) ~ 1, and imag{ris) ~ 0, the 
corresponding sturmians make a resonant contribution to the expansion. In Fig. [11] some of 
the rjs are illustrated in the form of an Argand diagram for three values of k, which shows 
that for s = 7 the values of r/y for all three values of k satisfy this resonance criterion by 
lying close to unity. 

III. APPROXIMATE QUASI PARTICLE METHODS Si AND ^2 

The original Q-P iteration procedure was developed for the solution of the L — S Eq. 
([3]), but is couched here in terms of a more general integral kernel O for a one dimensional 
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FIG. 9: Same as for Fig. d for s = 7 
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FIG. 10: Absolute value of ry^ as a function of the sturmian index s for two values of the wave 
number k. The eigenvalues rjg are defined in Eq. (jlOp with V = Vp. 
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FIG. 11: Argand diagram of rjs, for the potential Vp for s = 4, 5, ..10, for three values of k (in units 
of fm^^). For s = 7 they lie close to unity. 
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integral equation 

ilj{r) = F{r) + I 0{r,r')ij{r') dr'. (17) 





This equation is to be solved for ip, and F is a known "driving" term. The shorthand form 
of Eq. (dZ]) is 

ilj = F + Oip. (18) 

One procedure, convenient if O has a semi-separable nature, as is the case for Green's 
functions in configuration space, is to expand the unknown function into a set of Chebyshev 
polynomials, as is done in the lEM method |18|, and obtain an answer to high accuracy 
without iterations. If however the operator is very complex and of large norm, so that a Born 
series would not converge, then an alternative method based on expansions into sturmian 
functions, described by Weinberg [isl as the " Quasi-Particle" method, is as follows. One 
approximates the operator (9 by a separable part of rank A^, On, and the remainder Or is 
defined as 

Or = 0- On, (19) 

and iterates on the remainder. If the norm of On is less than unity, the iterations should 
converge. Since the numerical complexity of performing iterations is less than solving a 
linear equation with a matrix of large dimension, Weinberg's method can be computationally 
advantageous. Weinberg's method to obtain On is to calculate the sturmians Fj associated 
with the operator O which satisfy the same boundary conditions as ip 

or, = 7,r,. (20) 

and approximate (9 by a separable expansion into those sturmians Fj whose eigenvalues 7^, 
i = 1,2,.. ,N have a magnitude larger than unity. Iterations on the remainder will then 
converge, since they involve powers of 7^ which for i > N are less than unity. A numerical 
example (example P) will be given for the case that O = QV, where the potential V = Vp 
of Eq. (JT5II has a repulsive core and an attractive well. 

Since the computational effort to obtain the eigenvectors and eigenvalues of a general 
operator O can be larger than the effort to solve Eq. f|T7|) directly, a simplified method 
(method Si ) will be described below, that uses auxiliary sturmian functions. These methods 
avoid the rigorous calculation of the Fj by expanding the operator O into a finite set of 
auxiliary sturmian functions s = 1,2,.. .,N, which satisfy Eq. (JTOl) where V is an 
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auxiliary potential, and which obey the same boundary conditions as the function if). 

^ 1 

These sturmian functions can be obtained reliably and with high accuracy, using a spectral 
expansion method [l^, fl2|. However, the auxiliary sturmians are not used to obtain the 



eigenfunctions for the operator O, as is done in other studies ^\ where the emphasis was in 
obtaining the scattering resonances, but rather to set up a well converging iteration scheme. 
In method Si , designed to solve Eq. f ll7p iteratively one first defines a function J^*^^) which 
is the solution of 

J-f ) = F + 0^ J-f \ (22) 
Since On is given by the separable representation (12T]) . a solution of Eq. (!22|) is of the form 

TV 

^r = F + 5^ 0^s)c,. (23) 

s'=l 

with the coefficients c, defined as 



1 ,^ —(AT) 



s = l,2,..,iV. (24) 



These coefficients are obtained algebraically by solving the matrix equation (125|) for the Cg 
^ 1 

Y,{Ks'-Bs.s') c,= 3 = 1,2,.., AT. (25) 

s'=l ^ 

with 

This equation is obtained by multiplying Eq. fl2^ on the left by $s(r)\^(r)/ ($3^$^) and 
integrating over r from to 00. The iterative solution of Eq. ([1] 



V^ = J-r + x(^) + x(^^+x(^^ + ... (27) 

is obtained by solving 



X 



(n+l) 



+ n = l,2,... (28) 



iteratively with x*^^) = J-'}^''. This series will converge if the norm of On is less than unity. 
The advantage of this procedure is that the functions x^""* can also be obtained algebraically, 

13 



using the same matrix as required for the calculation of the coefficients Cg, and the evaluation 
of Orx^^^ only involves integrals over known functions, rather than the solution of linear 
equations, which would be more computer intensive. The matrices are of dimension N x N, 
where is of order 50 or less. Numerical examples are given below. 

In the approximate Q-P method iSg one solves the once iterated form of Eq. fllSp ip = 
F + 0{F+OiIj) 

ilj = F +OF+0^ij. (29) 

The iterative treatment is similar to that of method Si, with the driving function F replaced 
by F + OF, and the operator O in Eq. ( ITS]) replaced by O^. This method lends itself to 
a deconiposition of the matrices of in terms of the singular value decomposition method 
(SVD) [22|, [23|, that takes the place of finding the eigenvalues and eigenvectors 7^ and Fj 
of the operator O, and provides a measure of the norm of This method is described 



in Ref. 



24| . but is not used for the numerical calculations described in the text. 



IV. NUMERICAL EXAMPLES FOR METHODS Si AND ^2 

The operator 0{r,r') for all of the three numerical examples described below is 
QQ{r,r')V{r'), where Qq is the Green's function defined in Eqs. @ and (|5]), V is the scat- 
tering potential given by Vp, defined in Eq. ( IT5l) and Table [H The integral equation being 
solved is the conventional Lippmann-Schwinger (L — S) scattering equation, Eq. ([3]). For 
most cases the wave number k that defines the energy E = in Eq. ([2]) has the value 
k = 0.5 fm~^. As illustrated in Fig. [H Vp has a repulsive core near the origin followed by 
an attractive well that decreases exponentially with distance. The original Q — P method is 
illustrated by example P, and the approximate Q-P method S2 is illustrated for two sturmian 
potentials and Vws- The size N of the sturmian basis is 46, 50, and 41 for potentials V^, 
Vp, and Vwsi respectively 



A. Q-P example "P" 

In this case the sturmian potential is equal to the scattering potential, V = Vp, and 
hence the sturmian functions $3, defined by Eq. f lTOj) are the same, to within a normalization 
constant, as the eigenfunctions F of the operator O = QqVp, and the eigenvalues rjs are also 
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the same as the eigenvalues 7f'^^ The $s used in the present examples are normalized such 
that asymptotically they equal the function H{r), with unit coefficient. If an eigenfunction 
$s is normalized differently, then it can be renormalized according to 

where (FV^s) = /q'^™'''' -^(^)^('") '^^(r) dr. The value of k is 0.5/m~^, and r^ax = 30 fm. 
Some of these sturmians are illustrated in Figs. [6] to [9] and the absolute values of ris are 
displayed in Fig. [TOl The operator On becomes 

N 



=1 



In the present example = 26, and |?726| = 0.1425. The matrix elements in Eq. (125|) become 
diagonal, the coefficients c^, defined in Eq. fl25|) . are given by (Fy$s)/[(1 — rjs){^sV^s)]-i 
and the function J^(^) , solution of Eq. (122|) , is given by 



N 



The quantity 5*1 is defined by the asymptotic limit of J-" 



^r'w^n.)+Er^^*.«. (32) 



J^[^\r ^ oo) = F(r) + Si H{r). (33) 
and is given by 5*1 = , where 



and the iterative corrections to i/j are 



(34) 



^ = J^[''^+0rF+01F+0%F+ (35) 

as shown in the Appendix A. 

The values of q are shown in Fig. [121 They begin to decrease nearly monotonically 
for i > 15, but the rate of decrease is slow. For example cis = —0.7017 + 1.0067?, and 
C26 = —0.0107 + 0.0395?. This shows that without the additional iterations described next, 
the accuracy of S given by alone, would be good to only one or two significant figures. 
For example, the value of Si obtained from J^(^) is Si = 0.7075982522 + 0.493097515?, while 



the "exact" value of S, obtained via the spectral Integral Equation Method (lEM) [18| is 
S = 0.34016623828 + 0.8664518117?. In Fig. M the improvement due to the successive 
Weinberg iterations are displayed. 

15 



o 


• real 




o imag 


• 












• 





5 10 15 20 25 30 

sturmian index i 



FIG. 12: The coefficients Cj as defined in Eq. ()34p for tlie Q-P metfiod. Tlie sturmian potential is 
the same as the scattering potential Vp and k = 0.5 fm~^. 
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FIG. 13: The iterative improvement of the asymptotic value of S of the scattering wave function 
■0 is shown as the absolute value of the difference between the iterative value of S, obtained from 
Eq. (jSSp and the value of tp, obtained by solving the Lippmann-Schwinger scattering equation 
directly. The first point illustrates the difference between the asymptotic values of and while 
the subsequent points illustrate the contributions of the successive Weinberg iterations. 

B. Examples for the approximate Q-P Methods Si and 52- 

For both methods Si and ^2 the operator O is again QoV, with V = Vp, and the sturmian 
potential V can be either V^, Eqs. ( |T5l) . or Vp, the latter defined in Appendix B. In this 
section the normalization of the sturmian functions is such that 

{<^,V<^s)=Vs ,s = l,2,..,N (36) 

For method Si the elements of the matrix B, Eq. (126|1 . become Bs^s' = {^sV^s')Vs/ {^sV^s), 
and all other equations described in section III apply. The separable approximation to the 
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operator O, Eq. (12T]) . is given by 



0.-Y:S«Vt,)j^{<t,V. (37) 



and as a result of the normalization (136|) the n'th power of On is 



1 

(O^)" = J2 ^oV<!>s) {V "-^)^^, -($.'V^, (38) 

s,s'=l 

where the matrix V has the matrix elements 

V,,,, = ($.V^$.,). (39) 

Equation fl38l) shows that if the norm of the matrix V is larger than unity, iterations per- 
formed with {Om)"^ will not converge. For this reason separating from V a part that has a 
norm less than unity, and performing iterations for {OnY on that part, will converge. That 
is the reason for solving the iterated equation fl29l) rather than the original Eq. ([3]), as will 
be explained further below. For more complicated integral kernels a similar separation is 
also feasible with the Singular Value Decomposition method (SVD) [24] but is not needed 
for the examples given here. 

In method ^2 one decomposes O"^, the square of the full operator, into two parts O"^ = 
0% + 0\. According to Eq. (JSH]) {OnY is given by 

^ 1 

and one 0\ is defined as 

Ol = 0^- 0% (41) 

In order to solve the once iterated L — S equation (!29i) . one defines the function J^2 as 
the solution of 

J-f ^ = F + OF+Ol4''^ (42) 

given in the present case by 

AT 

J-f ) = F + GoV<^s)Vs,s'C^;^ (43) 

s,s'=l 
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where c^f) = {^s'VT2^^) /rjs'is the solution of the hnear equation 

J2(^s,s' - Vl,)c^;^ = -{<l>sV{F + OF)). (44) 

s'=l 

The terms X2"'* required for the subsequent iterations, 

^=^r+xr+xr+xr^+... (45) 

are obtained by solving 

X^2^'^ = Olxt'^ + Ol^-\ n = 1, 2, ... (46) 

with x^^) = ^ 

A justification for the advantage of method ^2 over S\ is as follows. Formally, Eq. f l46|) 
can be written as X2"^^'' = [(1 ~ ^at)^^^!;]" -^2^"*) ^^id hence the rate of convergence of the 
iterations depends on the norm of the operator (1 — 0\)~^0\. Since the norm of (1 — 0\j)~^ 
is expected to be smaller than unity, the rate of convergence should be faster than the powers 
of the norm of 0\. In view of the completeness of the sturmian functions, 0\ is also given 
by 

00 ^ 
0\= Q^y^s)V,^, — i,^,V. (47) 



By defining V 



V s,s> = Vs,s' s,s' > N 



and keeping in mind that 



Vs, s' = s,s' <N (48) 



— {^,,VgoV^s) = Vs',s (49) 

Vs' 



one can show that the powers of O]^ are given by 

00 ^ 

{Oir= E GoV^s)(v -i^s'V, (50) 

^ \ / S.S' Tla> 



s,s'=N+l 

and hence the norm of the iterations should decrease with the order of the iteration n as 



^2„-3 



n = l,2,.. (51) 



where t is some constant. 
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FIG. 14: Iterative rate of convergence of the asymptotic limit of the wave function to the "exact" 
one. The sturmian functions S and WS are obtained with potentials Vs and VwS: respectively, 
which have no repulsive core. Iteration method ^2 was used for results Vs and Vws and method Si 
was used for potential Vp. (The latter is identical to the potential used to calculate the scattering 
wave function). The solid line represents (0.66)^, with p = 2.6 x n — 3, n being the iteration order. 
The wave number is A; = 0.5 fm^^. 

By contrast, for method Si 

oo ^ 

{OnT= V ^o^$«)(f""') -{^s'V, (52) 

\ / s,s' Tl^i 

s,s'=N+l 

and by a reasoning similar to the one above, one expects 

n = l,2,.. (53) 

These convergence estimates are supported by Figs. [15] and [201 Although one would expect 
from Eq. f[5Tl) that the magnitude of the iteration error decreases as the square of the 
iteration order n, (i.e. a power 2) the fit in Fig. [15] decreases by a larger power, 2.6. 
Nevertheless, the convergence of method Si is considerably slower than for S2 , as illustrated 
in Fig. [201 and is in accordance with Eq. ( I53l) . 

The results for method 1S2, displayed in Fig. ( |T4l) . also show that the iterations using 
sturmians that are based on a potential with a longer range, Vwsi than the potential Vs 
converge faster. This feature can be understood from the fact that at large distances the 
sturmians become linearly dependent, since asymptotically they become proportional to 
each other, as can be seen from Figs. [1] to [7] and hence at large distances they do not form 
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FIG. 15: The convergence of the iterations n = 1, 2, .. as measured by the error of the asymptotic 
value of the wave function tp, Eq. ([6]). The results labeled B and B2 are obtained with sturmian 
potential V = Vb for methods S and 52, respectively. The solid line labeled as S2 is obtained with 
method ^2 for the sturmian potential Vs- The result P is obtained with the original Q-P method, 
with sturmians obtained for potential Vp. The open circles, labeled Greens, are obtained with 
Green's function iterations, based on Vp — Vb, described in Appendix B. 

a good basis set. Indeed, at short distances the convergence to the scattering wave function 
is much faster than at large distances, as can be seen from Fig. [16], obtained with sturmians 
for Vs. 

The results labeled B and B2 in Fig. [15] are obtained with methods S and 5*2, respec- 
tively, for sturmians based on potential Vb- This potential, defined in Eq. (fT6|) and further 
illustrated in the Appendix, is identical to potential Vp at large distances, but has its repul- 
sive core decreased by a factor that decreases to zero near the origin. Both potentials Vb 
and Vs decay in a similar fashion at large distances while Vb has a barrier near the origin, 
and Vs does not. The result that both sets of sturmians give nearly indistinguishable results 
for the iteration, as shown in Fig. [TS] shows that the behavior of the sturmians near the 
origin does not significantly affect the results. However, using the sturmians defined with 
the longer range potential Vws the convergence is significantly better at large distances, as 
is shown in Fig. [17] 

The open circles in Fig. [201 labeled as GreensB were obtained with a Green's function 
iteration method, described in Appendix B, in which potential Vp is divided into Vb + 
(Vp — Vp). The L — S equation with potential Vb is solved to produce the function 
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FIG. 16: Absolute value of the error of the wave function ip as a function of radial distance r . The 
result labeled "no iteration" illustrates the error of the function Eq. (gSD, with N = 46. The 

other line is obtained after the 14'th iteration. These resuts are obtained with method 52 using 
the sturmians for potential Vs- The scattering wave function ip is obtained with potential Vp for 
k = 0.5 fm~^. 
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FIG. 17: Same as Fig. [TBI using 41 sturmians based on the Vws potential. 



and the corrections due to {Vp — Vp) are obtained iteratively in a Born-series manner as 
approximation to the exact function ip. The asymptotic value of J-" is much closer to that of 
than for method S2, but the rate of convergence of the Green's function iterations is not 
as fast as that of method 
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V. SUMMARY AND CONCLUSIONS 



The main result of this paper is to present a viable iteration procedure to correct for 
the truncation errors of a sturmian expansion used for the solution of an integral equation. 
The basic idea is to construct a separable (and truncated) approximation to the integral 
kernel of the integral equation, that leads to an algebraic solution for the approximated 
integral equation, denoted as T . The corrections to T are performed iteratively, and the 
rate of convergence of the iteration is investigated. The advantage of doing iterations rather 
than solving the integral equation numerically directly, is that the numerical complexity for 
obtaining the algebraic solution for T and the subsequent iterations can be substantially less 
than for the direct numerical solution, especially if the kernel of the integral equation is very 
complicated. An example would be the solution of a Schrodinger equation with nonlocal 
potentials, recast into a Lippmann-Schwinger (L — S) integral form. The numerical examples 
are done in configuration space for a one-dimensional {L — S) equation for which the integral 
kernel is of the form QqV, where Qq is the undistorted Green's function and V the scattering 
potential, and the scattering energy is positive. However the formalism can be generalized 
to more general integral kernels. Several methods are compared. One is the "classical" 
quasi-particle 13| method, that requires a basis of sturmians that are eigenfunctions of the 
kernel of the integral equation. For the application of this method to the scattering of a 
particle from a potential with a repulsive core, Vp, the iterations converge very fast (denoted 
as method P) and to high accuracy, as illustrated in Figs. [131 CHI EDI If the sturmian 
eigenstates of the above mentioned integral kernel are not available, as is assumed here, then 
a set of auxiliary sturmian functions are introduced, and the convergence of the iterations 
is examined. Two methods are presented for this case. In method Si the iterative solution 
of the original integral equation is examined, while in method ^2 the solution of a once 
iterated integral equation is considered. It is predicted and also found that the iterations of 
method ^2 converge faster than for method Si , and both methods converge faster than the 
iterations performed with a Green's function Born approximation procedure, but more slowly 
than the original Q-P method. A key condition for obtaining an accurate approximation to 
the wave function at large distances is that the auxiliary sturmian functions be based on an 
auxiliary potential whose range is sufficiently larger than that of the scattering potential. 
An alternate method for solving the once iterated (but complicated) integral equation, based 
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on the singular value decomposition of a matrix, is described in Appendix C but not used 
in the numerical examples. 

An accuracy of at least eight significant figures is obtained for all the numerical calcula- 
tions performed. This accuracy is achieved by using the spectral method for solving integral 
equations 



19|, both for the scattering wave function and for the sturmian functions. As a 
result, previous investigations with positive energy sturmians |21|, hampered by lack of this 
type of accuracy, were improved upon and expanded in the present study. 

Appendix A 

A proof of Eq. fl35|l is as follows: The operator On, given by Eq. fj3T|) . projects any 
function into the space spanned by the sturmians s = 1,2,.., A^, while the operator Or 
projects any function into the sturmian space for s = A^ + 1, A^ + 2, ... The latter property 
follows from the orthogonality given by Eq. (|T2l) . and the fact that the sturmians form a 
complete set. Hence Or, Eq. (fT9|) . can also be written as 

oo 

Or =J2'^s[Vs/{'^sV<i>smsV (54) 

N+l 

This expression for Or also shows that the norm of Or is less than unity if |^7s| < 1 for 
s > A^. In view of Eqs. ([22]), and ([51]) one finds that O/jJ"}^^ = OrF. By writing 

^ = + X (55) 

and by by making use of Eqs. ([28]) for n = 1, (1- On)x^^^ = OrJ^^^^ = OrF. Hence x^^^ 
lies entirely in the sturmian space with s > N and x^^^ = OrF. In view of Eq. ( !28l) the same 
argument can be made for x^'^\ X^^'^ ■■ ^i-nd Eq. (15^ is proven for the case that O = GqV and 
V = V. 

Appendix B: The Green's function iterations. 

At large distances potential Vr is identical to Vp, but at short distances the repulsive 
core is reduced by the factor /(r) = 1 — exp[(r/0.5)^], as given by Eq. (fT6]) 



VB{r) = f{r) X Vp{r) (56) 

This potential is illustrated by the open circles in Fig. [18] This potential was defined for the 
purpose of performing a Green's function iteration, as follows. The potential Vp is divided 
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FIG. 18: The potental Vb is illustrated by the open circles, and is defined in Eq. ([56]) . The 
purpose of this potential is to eliminate the repulsive core of potential Vp, the latter illustrated by 
the + symbols. The difference between these two potentials, illustrated by the solid line, serves to 
perform iterations with a Green's function Born-type method. 

into two parts Vb and V2 = (Vp — Vb). The first part has a suppressed repulsive core, 
and the second part, V2 adds the repulsive core again. The iterations are based on Green's 
functions that explicitly include the distortion due to Vb, and iterate over the effect of V2 
in a fashion similar to the Born approximation. The Green's function method has many 
applications, but in view of the factor in the Greens functions F and H, Eq. (jl]), this 
method does not converge well {22! for low values of the wave number k. 




Appendix C: Iterative method to calculate Sturmian 
eigenfunctions. 



In what follows the eigenvalue subscripts s will be dropped, and the iterations which 
lead to the solution for a particular fixed value of s will be described. In the differential 
equation (|8]) the exact values of A and $ are initially unknown, but the equation can still be 
solved for a guessed value A*^^). One solutions, denoted as Y(r) satisfies the correct boundary 
conditions at the origin, and is integrated from inside outward. 



The other solution, denoted as Z{r), satisfies the correct boundary conditions asymptotically. 



and is integrated inward towards the intermediate point Rj chosen to lie within the region 
where V is non negligible. The three solutions at the point Rj are denoted as $/, Y/,and 




VY = A^^V Y, 0<r<Ri 



(57) 



V Z = A(i)V Z, Ri<r< 00, 



(58) 
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Zi. One can renormalize y by a factor t =Zi/Yi 

Y{r) = t Y{r) (59) 

so that Yj = Zj. 

If one multiplies Eq. (jH]) on both sides by Y, and Eq. ( 1571) by $, integrates over r from 
to Rj, and subtracts the results one from the other, after an integration by parts one obtains 

rRi 

(Y $' - Y' $) = (A - A^^^) / Y{r) V{r) $(r) dr, (60) 

Jo 

where the prime denotes a derivative with respect to r. By a similar procedure one obtains 

POO 

{Z $' - Z' $) Ig; = (A - A^i)) / Z{r) V{r) $(r) dr. (61) 

Jri 

By adding Eqs. fl60l) and fl6T]) . by noting that (F $' — $) |o = because Y and $ vanish 
at the origin, and that (Z $' — Z' $) |°° = because Z and $ obey the same boundary 
condition asymptotically (to within a normalization constant), one obtains 



$j(Z; - YD = (A - A(i)) 



F(r) V{r) $(r) (ir + / Z(r) l^(r) $(r) c?r 



(62) 



The above equation is still rigorously valid. The approximation consists in replacing $(r) in 
the first integral by Y{r), and by Z{r) in the second integral, and further, by replacing $7 
in the left hand side by either Yj or Zj. After making these approximations, and by dividing 
both sides of Eq. fl62|) by {YjY = {Zj)"^, one obtains the final result 

A(2) = aW + Z//Z1-Y//Y1 ^g^^ 

Y V{t) dr) /Y,^+ [j^^ Z \r) V{r) dr) /Z/ 

In the above, A*^^) replaces A as the iteratively corrected value for A^^^; also the factor t 
disappeared because it canceled in the numerator and denominator of Eq. ( l63ll . 

Equation ( 163|) is the generalization to sturmian eigenvalues of the iterative calculation of 
energy eigenvalues, Eq. (7) of Ref. [20]. It can also be applied to negative energies, and can 
be generalized to the cases of coupling potentials or of non-local potentials, which however 
goes beyond the scope of the present study. 

The functions Y and Z can be calculated as the solutions of the differential Eqs. (!57l) and 
( 158|) by any convenient finite difference method, or they can be obtained as the solutions of 
the integral equations 

fRi 

Y{r) = F{r) + A^^U go{r,r')V{r')Y{r')dr' (64) 
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and 



POO 

Z(r) = /7(r) + / go{r,r')V{r')Z{r')dr' (65) 

Jri 

or else the whole iteration procedure can be bypassed by obtaining the eigenvalues ris and 
eigenf unctions $s of the integral operator Qq V contained in Eq. fllOl) . In the numerical 
calculations described below, this latter procedure is carried out (with questionable accuracy 



according to Ref 



12|) so as to provide the initial values Ai^'' for the iteration, Eq. fl63|) . The 



present implementation obtains the solution of Eqs. and ( |65|) by using the spectral 



Chebyshev expansion method, [18|, 19| because its high accuracy can be predetermined by 
the specification of an accuracy parameter tol, and the method is by now well tested. 

Appendix D 

For future benchmark purposes, the eigenvalues of Eq. ([8]) are given in the Table [Tll The 
potential is defined in Eq. ( |T5i) and Table HJ the wave number is k = 0.5/m~^. 
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